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Abstract. Serious discrepancies have recently been observed between predictions of stellar evolution models in the 
0.7-1.1 Mq mass range and accurately measured properties of binary stars with components in this mass range. 
We study one of these objects, the eclipsing binary UV Piscium, which is particularly interesting because Popper 
(1997) derived age estimates for each component which differed by more than a factor of two. In an attempt 
to solve this significant discrepancy (a difference in age of 11 Gyr), we compute a large grid of stellar evolution 
models with the CESAM code for each component. By fixing the masses to their accurately determined values 
(relative error smaller than 1% for both stars), we consider a wide range of possible metallicities Z (0.01 to 0.05), 
and Helium content Y (0.25 to 0.34) uncorrelated to Z. In addition, the mixing length parameter qmlt is left as 
another free parameter. We obtain a best fit in the Teff-radius diagram for a common chemical composition {Z, 
y) = (0.012, 0.31), but a different MLT parameter aMLT,A= 0.95±0.12(statistical)-|-0.30(systematic) and QfMLT,B = 
0.65±0.07(stat)-f0.10(syst). The apparent age discrepancy found by Popper (1997) disappears with this solution, 
the components being coeval to within 1%. This suggests that fixing ttMLX to its solar value (~1.6), a common 
hypothesis assumed in most stellar evolutionary models, may not be correct. Secondly, since q:mlt is smaller 
for the less massive component, this suggests that the omlt parameter may decrease with stellar mass, showing 
yet another shortcoming of the mixing length theory to explain stellar convection. This trend needs further 
confirmation with other binary stars with accurate data. 

Key words. Stars: individual: UV Psc - Stars: fundamental parameters - Stars: eclipsing binaries - Stars: abun- 
dances - Stars: Hertzsprung-Russell (HR) diagram 



1. Introduction 

It is of the utmost importance for stellar evolution the- 
ory to match at least the best known objects (the Sun 
and non-interacting binary systems) before any attempt 
to derive proper information for star clusters or stellar 
populations in galaxies is made. 

Some binary stars are known to provide uncomparable 
astrophysical laboratories to calibrate theoretical stellar 
evolutionary models. Since the individual components of 
well-detached binary systems can be assumed to be two 
single-like stars with a common origin, they share the same 
chemical composition and same age, and, therefore, the 
observed parameters of both components are expected to 
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be matched by a single isochrone at the same chemical 
composition. 

Whilst in general a good overall agreement is found, 
some studies (pioneered by Popper (1997), see also e.g. 
Pols et al. (1997), Lastennet et al. (1999), Clausen et al. 
(1999) and Lastennet & Valls-Gabaud (2002) for a sys- 
tematic analysis) have unambiguously found that systems 
with components in the 0.7 to 1.1 Mq mass range raise 
serious difficulties. 

There are many well-studied systems with components 
in this mass range which present puzzling discrepancies 
with the predictions of stellar evolution theory. For ex- 
ample, if the seismic observations of a Cen A are fitted 
to state-of-the-art structure models, the resulting masses 
differ from the ones derived from dynamical analyses 
(Pourbaix et al. 2002, Thevenin et al. 2002). The visual 
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binary 85 Peg, some eclipsing binaries in the field, and 
in the Hyades open cluster are other clear examples (see 
Lastennet et al. (2002) for a brief review of the results 
obtained so far on these objects). 

This situation is very puzzling because one would ex- 
pect these stars to be reasonably well understood. Indeed, 
stars with masses larger than ~0.6 Mq bypass difficulties 
in the treatment of the equation of state and the atmo- 
sphere, while stars with masses larger than about 1.1 Mq 
have a permanent convective core, introducing an addi- 
tional parameter, the amount of overshooting, for their 
modelling. Hence, current stellar evolutionary models are 
expected to be able to match the basic properties of stars 
in the 0.7 - 1.1 Mq mass range, provided a good descrip- 
tion of the convection in their envelopes is used. 

One of these puzzling objects is the eclipsing binary 
UV Piscium (HD 7700, HIP 5980, hereafter UV Psc), 
a detached main-sequence binary (with components of 
types G and K), which is also a short period (<1 day) 
RS CVn^. This binary deserves particular attention ever 
since Popper (1997) derived age estimates for each com- 
ponent which differed by a factor larger than two, with an 
absolute age difference larger than 11 Gyr. Montesinos et 
al. (1988) classified its components as "apparently normal 
main-sequence stars" . What is the origin of this discrep- 
ancy? 

Popper (1997) used the Geneva models (Schaller et 
al. 1992) for Z =0.02 and Y =0.30, i.e. with a helium 
abundance to metallicity ratio fixed to AY / AZ=3, and 
aMLT fixed to a solar-calibrated value (q;mlt=1-63). 

Since the RS Canum Venaticorum stars show activ- 
ity, mass transfer episodes between active regions of both 
components are possible. In this case, the application of 
the single-star theoretical models used by Popper may not 
be relevant since the evolution of both components would 
not be independent any longer. However, for UV Psc, mass 
transfer does not account for the age difference because the 
computation of the Roche lobe radii (e.g. with the formu- 
lae given by Eggleton 1983) shows that none of these stars 
have overflown their Roche lobes. Therefore, we consider 
that UV Psc is a well-detached system and assume that 
the evolution of each component is independent. 

Another concern about RS CVn stars is the possible 
distortion of light curves due to starspots, which could 
bias the derivation of their fundamental stellar parame- 
ters. Jassur & Kermani (1994) suggested the presence of 
cool spots on the surface of the primary component, but 
according to Popper (1997, and references therein), this 
should not affect the value of the quantities used in this 
paper. The period has been recently revised by Sowell et 
al. (2001) who didn't find any evidence for a cyclic mod- 
ulation or period change. Hence activity cannot account 
for the age discrepancy either. 

^ RS CVn are binary systems where one of the components 
presents an extreme form of solar-like activity with starspots 
and variable magnetic fields. See Montesinos et al. (1988) for 
a discussion. 



To tackle this problem, we explore the influence of each 
physical parameter. A stellar evolutionary code allowing 
to compute tracks for different and independent physical 
parameters is needed for such a purpose, and so we choose 
the CESAM code (Morel 1997) which fulfills this condition 
and which has already been successfully tested in a variety 
of astrophysical contexts (e.g. Zahn et al. 1997, Suran et 
al. 2001, Cordier et al. 2002). 

In this paper, we present a detailed study of the system 
UV Psc with the CESAM models. Our goal is twofold: 1) 
to check if there exists a solution giving consistent ages for 
both components for a set of {Z, Y), and 2) to constrain 
aiviLT for both components and check if these values are 
consistent with the generally assumed solar value, aMLTo- 

The paper is organised as follows: Sect. 2 deals with 
the description of the CESAM models used and the grid 
computed for the purpose of this paper. Sect. 3 presents 
the results which fit both components of the eclipsing bi- 
nary UV Psc, with particular emphasis on the constraints 
derived from the Toft-radius diagram on the ajviLT param- 
eters. We also discuss the robustness of the results against 
changes in the Teff scale and in diffusion. Finally, Sect. 4 
provides a summary of our results. 



2. Evolutionary stellar models 

The stellar evolution calculations were computed with 
the CESAM code (Version 3, Morel ITMTjl . running at 
the Coimbra Observatory. Details on the physics of these 
models can be found in Lebreton et al. H1999|l . Here we 
just provide a short summary of the main CESAM phys- 
ical inputs: the CEFF equation of state is used, includ- 
ing Coulomb corrections to the pressure (Eggleton et al. 
119731 Christensen-Dalseaard [T^91|l : the nuclear reactions 
rates are from Caughlan and Fowler (1988); the solar 
mixture is from Grevesse and Noels (,1993) : the OPAL 
opacities (Iglesias and Rogers 199Sll used and com- 
plemented at low temperatures by opacity data from 
Alexander and Ferguson (1993) following a prescription 
by Houdek & Rogl (1996); the atmosphere is described 
with an Eddington r(r)-law; the convection is treated ac- 
cording to the mixing-length theory from Bohm-Vitense 
()1958|) . with the formalism of Henyey et al. (1965), giv- 
ing the mixing-length scale {aMurxHp) in terms of a free 
parameter times the local pressure scale height Hp. 

With these physical ingredients, we computed a solar 
model that fits the observed luminosity and radius with 
an accuracy of 10~*, with aMLT=l-63, helium abundance 
10= 0.268 and Zq~ 0.0175 for the commonly accepted 
solar age of ~4.6 Gyr (Dziembowski et al. [19913) and solar 
abundances as derived by Grevesse & Noels (|1993jl . 

Given this calibration, we computed a large grid 
of CESAM models from the Zero Age Main Sequence 
(ZAMS) to the terminal Age Main Sequence (TAMS, cen- 
tral hydrogen exhaustion), or for ages t younger than 10 
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Fig. 1. Location of UV Psc A and B (error bars in bold, the primary component being the hottest one) in the 
effective temperature-radius diagram with evolutionary tracks computed from the CESAM code for various chemical 
compositions. On each panel, tracks with the mass of the primary (dashed lines) and the secondary (solid lines) 
components are shown for 3 different values of omlt : Q:mlt=1-7, 1.2 and 0.8 (from left to right). For Z =0.01 and 
Z =0.02, we also show q:mlt=0.1 (solid line at the coolest temperatures) for the secondary component. All tracks are 
computed from the ZAMS to the TAMS or to 10 Gyr (see text for details). 



Gyr^. This was done for each component of UV Psc by fix- 
ing the masses to their accurately measured values: Ma = 
0.975 Mq and Mb = 0.760 Mq (Popper 1997). This ap- 
proximation is justified since the relative mass accuracy is 
better than 1% in each star. 

We considered different values of the metallicity Z (0.01, 
0.02, 0.03, 0.04 and 0.05), helium content Y (0.25, 0.28, 
0.31 and 0.34) and mixing length parameter ofMLT (0.8, 1.3 
and 1.7) in a large three-dimensional grid, since we do not 
want to introduce correlations between these parameters. 



UV Psc being a binary belonging to the close solar neigh- 
bourhood (d ~ 63 pc), ages older than the Galactic disc (likely 
to be 8-10 Gyr old, at most, see Carraro 2000) should be irrel- 
evant here. 



We also computed additional models for the secondary 
component with aMLT=0.1 for Z = 0.01 and Z = 0.02^ 
We stress that the quantities Z and Y are independent, 
unlike most of the extant publicly available theoretical 
models where Y is derived from a fixed law once Z is de- 
termined. All the models computed with the input physics 
described above are displayed on Fig. 



^ The internal structure of such Iow-cimlt models does not 
present any convective region at all (this can be explained by 
the fact that the convective efficiency is proportional to the 
square of the mixing length and hence to o^lt)- The discussion 
on how realistic these models are is beyond the scope of this 
paper, but we would like to point out that we only used them 
to interpolate models with omlt slightly lower than 0.8. 
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3. Analysis of the results 

UV Psc is an eclipsing binary with very accurate masses 
and radii derived for both components. The fundamental 
data used in this paper (from Popper, 1997) are listed in 
Table 1 and are a key point of this study: the measures by 
Popper represent such a significant revision and improve- 
ment that a comparison with older data is irrelevant. 

For instance, Montesinos et al. (1988) provided a com- 
prehensive work on RS CVn EBs and derived a metal-rich 
solution for UV Psc {Z =0.04 and Y =0.25) from the log 
g-Teff diagram, but the data they used were very signifi- 
cantly revised by Popper (AM^=0.3 M©, etc) hampering 
any useful comparison. 



log R{track) — log Ri 
^(bii?0 



(1) 



We assume that both components share the same chemical 
composition, i.e. Za = Zb and Ya = Yb- Moreover, they 
should also have equal (or at least consistent) ages, hence 
\\tA — is II should be minimal. 

In Fig. 13 we show all the solutions in three projec- 
tions: Z-Y, Z— aMLT,A and Z—aMiJT,B diagrams. The 
solutions (small dots on Fig. which have a good for 
both components (xf < Xi,min + 1 ^'^th i =A, B) show 
that i) high metallicities {Z >0.032) are indeed ruled out, 
ii) 1.7 >aMLT,A > 0.9, and iii) 1.0 >aMLT,B > 0.6. 



Table 1. Parameters from Popper (1997) for the primary 
(A) and secondary (B) components of UV Psc. 



Comp. 


M 


R 


log (Teft/K) 




[Mo] 


[R©] 




A 


0.975±0.009 


1.11±0.02 


3.762±0.007 


B 


0.760±0.005 


0.83±0.03 


3.677±0.007 



First, before exploring all possible solutions in the {Z, 
Y, q;mlt, t) parameter space, interesting qualitative con- 
clusions can be derived from a simple inspection of Fig. 1. 
Whatever the values of Y and auur are, all models with 
Z>0.04 are unable to match simultaneously the location 
of both components, strongly suggesting that metal-rich 
solutions can be ruled out. In addition, while a fit to 
both components seems to be possible in some panels - 
e.g., the upper left panel at {Z,Y) =(0.01,0.25)- incon- 
sistent ages appear, excluding these possible solutions. 
This is also the case for the solution found by Popper 
at {Z,Y) =(0.02,0.30) using the Geneva models. 

Altogether, this means that we expect to find solu- 
tions with consistent ages for chemical compositions in 
the range Z =0.01-0.02 and Y =0.28-0.34. At this level of 
qualitative discussion, it appears difficult to infer strong 
conclusions about the omlt parameter but the panels of 
Fig. 1 suggest that aMLT,B may be systematically smaller 
that aMLT,A- 

To obtain more detailed conclusions, a systematic work 
matching some selected relevant criteria has to be per- 
formed. For this purpose, we explore by interpolation"* the 
{Z, Y, QfMLT, t) parameter space described in the previ- 
ous section to find out the best mathematical solutions 
matching the location of UV Psc in the Tcff-radius dia- 
gram for both components (i = A, B). This is performed 
by minimizing the x^-functional 



X^{Zi,Yi,aMUTi,ti) = 



logT,s{track)-\ogT, 



effi 



o- (log TcSi ) 

Z, Y and qmlt in steps of 0.001, 0.005 and 0.05, respec- 
tively. 
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Fig. 2. (Z, Y, aMLT,A, "mlt.b) solutions (dots) matching 
the location of both UV Psc A and B in the Tcff-radius 
diagram (ToffS from Popper, 1997) with the CESAM evo- 
lutionary tracks. Adding a constraint to obtain consistent 
ages for both components (At/t<20%) provides the sub- 
sample shown as crosses. The best agreement in age (0.9%) 
is shown with a filled circle. The Sun (dotted circle) is also 
shown for comparison. 



Adding the constraint that the ages are consistent 
with each other (At/t<20%) provides the subsample 
shown with crosses. It defines a much smaller region 
(Fig. 121) and shows that i) the metallicity is solar or 
sub-solar {Z <0.019) and the helium content is larger 
than the solar value (Y >0.28), u) 1.1 > aMLT,A > 0.9, 
and iii) 0.7 > ajviLT.B > 0.6. 

The best solution (filled circle in Fig. ^ is obtained 
with a Z =0.012, Y =0.31 chemical composition and pro- 
vides an excellent agreement in age (0.9% discrepancy). 
We predict an age of about 1.9 Gyr for both components 
of UV Psc A (1.92i°:^ Gyr) and B {IMtH Gyr). While a 
unique {Z, Y , t) combination is needed to simultaneously 
match the (masses, radii, Teffs) of both components, this 
solution implies a different value for the convection pa- 
rameter for the two components of the binary: omlt (UV 
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Fig. 3. Best fit solution (filled circle in Fig. 2) for both 
UV Psc A and B in the temperature-radius diagram (Toff s 
and radii from Popper, 1997) with the evolutionary tracks 
computed from the CESAM code: aMLT,A=0.95 ±0.12 
and aMLT b=0.65 ±0.07. 



Psc A) = 0.95±0.12 and aMLT(UV Psc B) = 0.65±0.07. 
The fit provided by this solution is shown in Fig. 13 

The non-solar values for qmut parameters that we de- 
rive obviously deserve further comments, but the robust- 
ness of these solutions needs to be discussed first. 



3.1. Influence of the T^s scaling 

Since our main concern was to solve the age discrepancy 
found by Popper (1997), we used the same values for 
the mass, radius and T^g. We have shown in the previ- 
ous section that this discrepancy can be solved within the 
same data set. However, whilst the individual masses and 
radii are very accurately defined, the TeffS are not direct 
quantities, and depend on colour indices and photometric 
calibrations. We would like to point out that according 
Popper l|TM5|l the presence of dark spots has relatively lit- 
tle effect on the color index of a star. To test the possible 
influence of the T^g scaling on our results, we used as an al- 
ternative the metallicity-dependent empirical calibrations 
of Alonso et al. ijTMSIl . T^s =f(S-F, [Fe/H]) and used the 
B — V indices of both components from Table 6 of Popper 
(1997) assuming that the reddening can be neglected be- 
cause UV Psc is in the inner solar neighbourhood (^ 63 pc 
according to its Hipparcos parallax) . Since there is no di- 
rect determination of the iron abundance [Fe/H] available 
for UV Psc, we have computed both Tcff s from the Alonso 
et al.'s formula (1) assuming three different metallicities: 
Teff,A= 5574 K and T^s.b^ 4506 K (Z = 0.01), Teff^A^ 
5682 K and T^s^^ 4586 K (Z = 0.02), Teff' a= 5754 K 
and rcff,B= 4638 K (Z = 0.03). The uncertainty on the 
primary (secondary) component is less than 150 K (100 
K), and so - to be conservative - we consider these upper 
limits as the uncertainty. For further comparisons, we take 
the TeffS derived for Z =0.02. They are cooler by about 



100 K (170 K) for the primary (secondary) component in 
comparison to the T^ss that we have used previously (see 
Table 1). 

The influence of adopting the Tegs we have derived 
from the Alonso et al. calibrations can be seen in Fig. 01 
A comparison with Fig. |21 shows that the global solutions 
(dots) are fully consistent with each other. However, the 
best solution according to the criteria defined in Section 
3 is slightly different: (Z, Y) = (0.015, 0.295), the ages 
predicted being about twice older (UV Psc A: ~ 4.07 Gyr 
and B: ~ 4.14 Gyr). In spite of these shifts, we note that 
a different value is still needed for the convection param- 
eter of the two components (aMLT,A= 1-15 and aMLT,B= 
0.65). The difference AaMLT is even increased in com- 
parison to the results inferred from the Popper (1997) 's 
Teff scale, and hence the conclusions of the previous sub- 
section remain qualitatively unchanged. Quantitatively, a 
systematic error of about +0.20 in aMLT, A (and a neghgi- 
ble one in aMLT.e) may be due to the uncertainties in the 
temperature scaling. 
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Fig. 4. Same as Fig.|5]but with the Tegs derived from the 
Alonso et al. (1996) calibrations. The best agreement in 
age (1.6%) is shown as a filled circle. 



3.2. Influence of diffusion 

Another potential source of systematic error is the in- 
fluence of diffusion. We have treated this by comput- 
ing CESAM models incorporating the microscopic diffu- 
sion as described by the simplified formalism proposed 
by Michaud & Proffitt (1993) with metal elements as 
trace elements. Note that the radiative acceleration can 
be neglected for sub-solar stars (Turcotte et al., 1998). 
As shown on Fig. |31 diffusion almost mimics an increase 
in aMLT, even though the slopes of the tracks are obvi- 
ously changed. At fixed radius, the effect of diffusion is to 
change the effective temperatures by some 65 K (55 K), 
well within the l-a uncertainties in TeS discussed above. 
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Alternatively, taking these uncertainties at face value, one 
would have to increase aMLxby about 0.11 (0.07) at most. 

In summary, even though the absolute values of 
aMLT may change due to these combined effects, the trend 
of increased ajviLxwith mass remains robust. 
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Fig. 5. The effect of including microscopic diffusion mod- 
els, assuming the same composition (Z,Y)=(Q. 012, 0.31) 
for each component. 



More results are needed before reaching any definitive 
conclusion, because a larger sample may well reveal a scat- 
ter of c^MLT for a given mass. For instance, recent calibra- 
tions of q;mlt in visual binary stars were performed by P. 
Morel and collaborators using the CESAM code (with dif- 
ferent physical ingredients than those used in this paper) . 
For the stars in the range of mass relevant in this paper, 
they find a Cen B, M = 0.97 M©, q;mlt= 1-53 (Morel 
et al 2000a), i Peg B, M =0.81 Mg, aMLT=l-36 (Morel 
et al 2000b) and 85 Peg A, M =0.84 Mg, aMLT=1.80 
(Fernandes et al. 2002). Since both t Peg B and 85 Peg A 
have very similar masses, this seems to support the idea of 
a dispersion of omlt at a given mass. Solar models with 
and without diffusion yield different values in the inter- 
nal convective regions (1.7 vs 1.9, Richard et al. 1996). 
There are no further data for main sequence stars with 
larger masses. Studies of more massive stars have dealt 
exclusively with red giant stars. For instance, aMLT seems 
to increase to values up to 3 to reproduce the red giant 
branches in open clusters (Stothers & Chin, 1997), while 
Iwamoto & Saio (1999) present a trend of ofMLT increas- 
ing with mass and metallicity again for giant stars. In 
these stars overshooting is playing a major role and so 
disentangling the relative effects of mass, metallicity and 
overshooting is far from trivial. 

Therefore, at this stage, a relation between mass and 
ttMLT along the main sequence remains highly premature 
and must be seen with caution. 



3.3. Universality of the omlt parameter 7 

We have quantitatively confirmed the suggestion made by 
Clausen et al. (1999) that the adoption of a significantly 
smaller Omlt for the secondary component may remove 
(or at least decrease) the general problem of age discrep- 
ancy presented in the Introduction. While it is obvious 
that only observational constraints on the chemical com- 
position of UV Psc would allow a proper test, it seems 
that this problem can be solved theoretically for this bi- 
nary. The solution obtained for this system suggests that i) 
the components have a different Q!mlt (Aq!mlt=0.3±0.14) 
and ii) these values are clearly different from the solar one. 

This implies that the solar value of the Omlt is not 
universal, an hypothesis made in most of the currently 
available stellar evolutionary models. 

As for the trend with mass, a recent study on the 
Hyades eclipsing binary vB22 (Lebreton et al. 2001) seems 
to support it. We show in Fig. that UV Psc and vB22 
are clearly not consistent with the solar aMLT value. 
Moreover, Lebreton et al. (2001) show that the slope of the 
main-sequence in the Hyades suggests that aMLT could de- 
crease with mass, in particular they obtained aMLT<l-4 
below 1 Mq. The results we have obtained for UV Psc 
give additional support to this trend. However we note 
that this trend seems to contradict the results of Ludwig 
et al. (1999) coming from detailed hydrodynamical sim- 
ulations where, for a fixed gravity, aMLT increases with 
deacreasing Toff. 
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Fig. 6. Best fit solution for both UV Psc A and B in a 
mass-aMLT diagram. The value of omlt.a is sligthly dif- 
ferent according to the assumed TcffS: aMLT=0.95 (Popper 
1997), aMLT=l-15 (Alonso et al. 1996, see text). The so- 
lution is independent of the choice of the Tcgs for the 
secondary component. The Sun (cf. §2) is shown for com- 
parison, as well as estimates (shown as upper and lower 
limits) derived from CESAM models with the same input 
physics for both components of the Hyades binary vB 22 
(Lebreton et al. 2001). 



We also note, in passing, that our results are not in 
contradiction with the recent analysis by Palmieri et al. 
(2002), who on the basis of the colours of red giants in 
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globular clusters concluded that there is no dependence 
of aMLT with metallicity. The extremely narrow range in 
mass probed by these red giants would be nearly equiv- 
alent to a single point on our diagram. Since the de- 
tailed modelling of red giants is more complicated than 
the one of main sequence stars, a larger sample of double- 
lined, eclipsing binaries with components in the 0.7-1.1 
Mq would offer a much cleaner and definitive test. 

In any case, all those possible variations and their 
trends reflect a shortcoming of the MLT theory more than 
a real physical effect. If one insists in using MLT, its pa- 
rameter Q!MLT should be varied. This is perhaps the first 
time that such a result is obtained using observations, 
since the variablity of the effective aMLT is a well-know 
result of detailed hydrodynamical simulations of stellar 
convection (Ludwig et al. 1999). 

On the other hand we note that the value of omlt 
determined in atmosphere of the stars gives in general dif- 
ferent results in relation to those obtained from the stel- 
lar interior analysis. For instance, van't Veer-Menneret & 
Megessier (1996) found that the spectral fitting of the first 
for Balmer lines of the solar spectrum required aMLT= 0.5, 
a result confirmed by Fuhrmann et al. (1993) for other cool 
stars, while Gardiner, Kupka & Smalley (1999) found that 
Q!MLT ~ 0.5 for TeffS below 6000 K but aMLT > 1-25 in 
the range 6000-7000 K to possibly decrease again to omlt 
= 0.5 above 7000 K. 

3.4. Implications for the helium to metal ratio 

Since we have derived Z and Y independently for UV 
Psc, the helium to metal ratio (AY/AZ) of this system 
can be inferred assuming a primordial helium abundance. 
The initial helium content is fixed by taking the currently 
assumed law Y = Yp + (AY/AZ) x Z, where Yp is the 
primordial helium abundance. If we consider Yp =0.235 
(Peimbert et al. 2000) and use the {Z, Y) solution that 
we have found for the UV Psc binary, we derive AY/AZ ~ 
6 and 4 by using the Popper (1997) or Alonso et al. (1996) 
Tcirs, respectively. In both cases, the ratio is large in com- 
parison to what can be obtained from the solar model 
described previously: (AY / AZ)q ~ 1.9. This may sug- 
gest that the currently assumed simple linear relationship 
between Z and Y is not universal and/or that AY/AZ is 
not unique. 

The possible relation between Y and Z is one of the most 
interesting and puzzling problems in astrophysics since the 
pioneering work by Faulkner (1967). 

The most recent observations in HII regions and 
planetary nebulae as well as semi-empical determinations 
using theoretical stellar models - including the detailed 
calibration of the Sim - (cf. Ribas et al. 2000, Fernandes 
2001 and references therein) support values around 2 or 
3 , assuming Yp ~ 0.23. Nevertheless, some exceptions 
indicate values larger than 4 (e.g. Pagel et al. 1992 from 
observations of extragalactic HII regions, and Belikov 
et al. 1998 from the main sequence of the Pleiades). 



Our determination clearly supports the idea of a large 
AY/AZ ratio. 



4. Conclusions 

We study the eclipsing binary UV Psc, which is of partic- 
ular interest since Popper (1997) derived age estimates 
for both components which differed by some 11 Gyr. 
Adopting his same data set (masses, radii and Tegs) we 
compute a large grid of CESAM models for each compo- 
nent^. We obtain a best fit in the Teff-radius diagram for 
a common chemical composition {Z, y)=(0.012, 0.31), a 
negligible age difference, and a different MLT parameter 
(Q!mlt,a= 0.95±0.12 and aMLT,B= 0.65±0.07). The dif- 
ference aMLT, A — Q;mlt,b = 0.3 is significant and remains 
robust against changes in temperature calibrations and 
the inclusion of diffusion. 

Within the framework of the MLT theory, this suggests 
that fixing aMLT to its solar value (~1.6), the usual hy- 
pothesis made by most stellar evolutionary models, may 
not be correct. Secondly, since the aMLT is smaller for the 
less massive component, this suggests that the aMLT pa- 
rameter may decrease with the stellar mass. This trend 
needs further confirmation with other binary stars with 
accurate data. Again, this may refiect a problem of the 
application of the MLT theory, more than a real physi- 
cal effect, but it has to be taken into account when using 
MLT. On the other hand as the main observational errors 
come from the Tcff, new and independent determinations 
are needed for both stars. 

So far, few other binaries have data with the precision 
required to confirm this trend and observational programs 
are highly needed (e.g. Kurpinska-Winiarska & Oblak, 
2000). 
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